Creating Single Collective Atomic Excitations via Spontaneous Raman Emission in 
Inhomogeneously Broadened Systems : Beyond the Adiabatic Approximation 
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The creation of single collective excitations in atomic ensembles via spontaneous Raman emission 
plays a key role in several quantum communication protocols, starting with the seminal DLCZ 
protocol [L.-M.Duan, M.D. Lukin, J.I. Cirac, and P. Zoller, Nature 414, 413 (2001).] This process is 
usually analyzed theoretically under the assumptions that the write laser pulse inducing the Raman 
transition is far off-resonance, and that the atomic ensemble is only homogeneously broadened. 
Here we study the impact of near-resonance excitation for inhomogeneously broadened ensembles 
on the collective character of the created atomic excitation. Our results are particularly relevant for 
experiments with hot atomic gases and for potential future solid-state implementations. 
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I. INTRODUCTION 

The creation of collective atomic excitations via spon- 
taneous Raman emission plays an important role in well- 
known quantum information protocols, in particular the 
quantum repeater [l[ protocol proposed by Duan, Lukin, 
Cirac and Zoller (DLCZ) B- for other protocols based 
on the same process see Ufl IE HI- The basic (ideal- 
ized) scheme is as follows, cf. Fig. [TJ In an ensemble of 
three-level systems with two ground states g and s and 
an excited state e all N atoms are initially in the state g. 
An off-resonant laser pulse on the g — e transition (the 
write pulse) leads to the spontaneous emission of a Ra- 
man photon on the e — s transition (the Stokes photon). 
Detection of this photon in the far field, such that no 
information is revealed about which atom it came from, 
creates an atomic state that is a coherent superposition 
of all the possible terms with N — 1 atoms in g and one 
atom in s, in the simplest case the completely symmetric 
state 
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FIG. 1: Basic level scheme for the creation of collective atomic 
excitations in atomic ensembles via spontaneous Raman emis- 
sion. All atoms start out in g. A laser pulse off-resonantly 
excites the g — e transition, making it possible for a photon 
to be emitted on the e — s transition (with small probability) . 



Such a state corresponds to a single collective atomic ex- 
citation in s. In general the term with the n-th atom 
in s will have a phase e I ( k ™- k s)x„ j w here \t w is the k 
vector of the write laser, kg is the k vector of the de- 
tected Stokes photon, and x„ is the position of the n-th 
atom. The phases in Eq. |T]) thus correspond to the case 
ks = k w . Moreover in practice the amplitudes of the 
different terms may vary, depending on the laser profile 
and the shape of the atomic ensemble. 

A remarkable feature of such collective excitations that 
is of great interest for practical applications is that they 
can be read out very efficiently by converting them into 
single photons that propagate in a well-defined direction, 
thanks to collective interference @, 0, Q . Resonant laser 
excitation of such a state on the s — e transition (the read 
laser pulse) leads to an analogous state with N — 1 atoms 
in g and one delocalized excitation in e. All the terms 
in this state can decay to the initial state \g)® N while 
emitting a photon on the e — g transition (the Anti-Stokes 
photon). If the phase matching condition ks + ^-as = 
k w + k r is fulfilled, where k r is the k vector of the read 
laser and ItAs that of the Anti-Stokes photon, then the 
amplitudes corresponding to the various terms interfere 
constructively (provided that there are no other effects 
disturbing the interference, such as atomic motion or the 
effects studied in this paper, cf. below), leading to a 
very large probability amplitude for emission of the Anti- 
Stokes photon in the direction given by k. w +k r — ks. For 
atomic ensembles that contain sufficiently many atoms, 
emission in this one direction can completely dominate 
all other directions. This allows a very efficient collection 
of the Anti-Stokes photon @, Q. 

Note that there is no such collective interference ef- 
fect for the emission of the Stokes photon, since its emis- 
sion by different atoms corresponds to orthogonal final 
states, e.g. the state |s)|<?)...|<7) if the Stokes photon was 
emitted by the first atom etc. Full "which-way" infor- 
mation about the origin of the photon is thus stored in 
the atomic ensemble, making interference impossible [9j. 
As a consequence the total emission probability for the 
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Stokes photon is simply given by the sum of the emission 
probabilities for each atom, and there is no preferred di- 
rection of emission. 

The creation of collective excitations via spontaneous 
Raman emission is usually analyzed theoretically under 
the assumption that the write pulse is far off-resonance 
HBH ( but see Section VLB of Ref. [11,]) . Under this 
condition it is possible to adiabatically eliminate the ex- 
cited state. However, in experiments the far off-resonance 
condition is frequently not fulfilled. From an experimen- 
tal point of view it can be advantageous to approach res- 
onance in order to increase the rate for the spontaneous 
Raman process or in order to avoid exciting nearby lev- 
els. For excitation relatively close to resonance it is no 
longer justified to eliminate the excited state. The pre- 
cise frequ ency of the excited state then influences the 
dynamics [12j , and it becomes important to consider the 
effects of inhomogeneous broadening of the transition 
between ground and excited state [1 31 ] - which is signif- 
icant in many experimental situations, e.g. for hot gases 
[II EE El, El, where the relevant mechanism is Doppler 
broadening. (Note that the effects of Doppler broaden- 
ing are negligible in similar experiments with cold atomic 
gases 0,11 El.) 

Inhomogeneous broadening will also be an essential 
factor in future experiments with solid-state atomic en- 
sembles, in particular rare-earth doped crystals 
where it is due to the crystal environment. These systems 
are otherwise very attractive candidates for realizing the 
DLCZ and similar protocols thanks to their excellent co- 
herence properties. For example, storage times exceed- 
ing 1 second have already been demonstrated in such a 
system for coherent atomic excitations in s created via 
electromagnetically induced transparency [20] , and light 
at the single-photon level has been stored and re-emitted 
using the "atomic frequency comb" protocol [2l|,|22j. Ref. 
[l9| proposed to reduce the inhomogeneity in these solid- 
state systems via spectral tailoring techniques similar to 
those employed in light storage experiments [2fJ, |2jj, |23| • 
However, such an approach greatly reduces the number 
of available atoms, making it much harder to write and 
read the atomic excitations efficiently. 

Motivated by these considerations, we here analyze 
the creation of collective atomic excitations in inhomoge- 
neous systems by spontaneous Raman emission without 
resorting to the usual adiabatic elimination of the excited 
state. This makes it possible to quantify the impact of 
near-resonance excitation in combination with inhomoge- 
neous broadening on the collectivity of the created atomic 
excitation. We introduce the term collectivity for the fi- 
delity of the created excitation with respect to the ideal 
state of Eq. ^ . This quantifies the degree of collective 
interference that is possible when reading out a given 
atomic excitation. It is equal to one (corresponding to 
the possibility of perfect collective interference) for exci- 
tations that are created under far off- resonant conditions. 
It is reduced for near-resonant excitation in inhomoge- 
neous systems. In such systems, atoms closer to reso- 
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FIG. 2: Evolution of the absolute value of the excited state 
amplitude for a single atom (that starts out in the ground 
state) under pulsed excitation for different values of detuning 
A. The fine line shows the exciting laser pulse Q(t), which is 
a Gaussian centered at 0.2 /is with FWHM pulse duration of 
0.1 /Us and a maximum Rabi frequency of 1 MHz. The same 
pulse is used in the examples throughout the paper. 

nance with the write laser will have a larger amplitude of 
emitting a Raman photon. Moreover for near-resonant 
(non-adiabatic) excitation the excited state plays a role 
in the Raman process, leading to phases that differ from 
atom to atom in inhomogeneous systems. These effects 
perturb the collective interference that is at the heart of 
the read-out process. 

The main goal of our present work is the characteriza- 
tion of the collective atomic excitation that is created by 
the emission and detection of the Stokes photon, with a 
focus on spectral aspects due to near-resonant write ex- 
citation and inhomogeneous broadening. Spatial effects 
have previously been analyzed in detail in Ref. [101 ] , and 
propagation effects for the Anti-Stokes photon in Ref. 
[111 ]. This paper is organized as follows. In section II 
we study spontaneous Raman emission for a single atom 
under pulsed excitation. We show that the amplitude 
for detecting a Stokes photon at a given time is propor- 
tional to the amplitude of the atom being in the excited 
state. In section III we show how the single-atom results 
can be used to quantify the collectivity for an inhomo- 
geneously broadened atomic ensemble. In section IV we 
give numerical examples relevant to hot atomic gases and 
solid-state systems. Section V contains our conclusions 
and an outlook towards future work. 



II. SPONTANEOUS RAMAN EMISSION FOR A 
SINGLE ATOM UNDER PULSED EXCITATION 

Let us first recall the well-known dynamics of a two- 
level system with levels e and g under pulsed excitation, 
cf. Fig. [2] Suppose that the system starts out in g and 
that it is excited by a pulse with Rabi frequency O(i) 
that is detuned from resonance by A. We focus on the 
case J dtn(t) < 1 (moderate pulse area). Depending on 
the value of A there are different regimes. For far off- 
resonance excitation, A 3> the amplitude to be in 
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e is well approximated by — such that the system 
returns to g after the pulse. This is called the adiabatic 
regime, where the system always stays in the momen- 
tary lowest-energy eigenstate of the Hamiltonian, cf. the 
lowest curve in Fig. [2] For smaller values of A the popu- 
lation in e attains a maximum during the pulse and then 
declines, but without returning exactly to zero, cf. the 
intermediate curve in Fig. O Finally, for resonant ex- 
citation the population in e increases monotonously for 
the whole duration (as long as one can neglect the spon- 
taneous decay of the excited level), cf. the top curve in 

Fig.rj 

Now add a third level s. We are interested in sponta- 
neous emission on the e — s transition. This means that 
we have to include a continuum of vacuum modes of the 
electro-magnetic field. We follow the Wigner-Weisskopf 
treatment of spontaneous emission |24j. The relevant 
states are now |<?)|0), |e)|0), where |0) is the vacuum state 
of the field, and the set |s)|l w ), where \l u ) are a contin- 
uum of one-photon states labeled by the frequency u. 
For simplicity, and motivated by typical experiments, we 
focus on one specific direction of emission, treating the 
effect of the other modes globally as a contribution to 
the decay of the excited level e; note that e can moreover 
also decay on the e — g transition. 

Denoting the amplitudes of the states |g)|0), |e)|0) and 
|s)|l w ) by a, /3 and 7^ respectively, and the overall decay 
rate of e by T, the dynamical equations are 

a = -in/3 (2) 
f} = -ma-i&P-TP-igY^lu (3) 

iu, = -ig/3 ~ iwjuj- (4) 

The coupling constant g depends on the dipole moment of 
the e — s transition, but also on the quantization volume 
and on the solid angle of the mode under consideration. 
Its precise form is not important for our purposes here, 
cf. Ref. [H. The solution of Eq. g]) is 

lUt) = ~ig [ dtfe-W-Vptf). (5) 
Jo 

Inserting this into Eq. (3) gives 
$= -ina-iAf3-Tf3-g 2 V f dt'e'^-^ j3(t'). (6) 

Changing the order of summation and integration in the 
last term, and making a Wigner-Weisskopf type approx- 
imation, 

£V*"(t-f) ocJ(t-t'), (7) 

we see that the modes under consideration just make a 
contribution to the overall decay term — T/3, as was to be 
expected. The relative size of this contribution depends 
on the solid angle of the considered mode (and also on the 



branching ratio between the e—g and e—s transitions). It 
is essentially negligible in typical experimental situations 
where the solid angle of collection is small. More inter- 
esting for our purposes is the effect of the detection of 
a Stokes photon at time t. Since the annihilation opera- 
tor for the mode under consideration satisfies a — a u 
(in the Schrodinger picture), the amplitude for such a 
detection is 

c (t) = y>u,(t) = - y> f dt'e-^-^m « m, 

where the equality follows from Eq. (J5]) and the propor- 
tionality from Eq. ([7]) . This shows that the amplitude for 
detecting a Stokes photon is simply proportional to the 
amplitude of the atom being in the excited state. Below 
we will show that, as a consequence, the (numerical) solu- 
tion of the two-level problem gives us all the information 
we need in order to study the collective interference. For 
a single atom, the atomic state conditional on detecting 
a Stokes photon is simply \s) in the single-atom case. It 
is more interesting in the iV-atom case below. 

III. COLLECTIVITY IN INHOMOGENEOUSLY 
BROADENED SYSTEMS 

We now consider the situation where the laser pulse 
excites an ensemble of N atoms, which do not all have 
the same resonance frequency. Without the third level s 
one would have a time-dependent state 

N N 

H(a n (t)\g) n + f3 n (t)\e) n ) = JJ |G n (t)>. (9) 

n— 1 n— 1 

We are interested in the case where there is a third level, 
but where the spontaneous emission of a photon on the 
e — s transition into the considered directional mode oc- 
curs with only a small probability. We then take into 
account only terms that correspond to a single emission. 
Detection of a single Stokes photon at time t creates a 
conditional state proportional to 

Cl (t)\ s )\G 2 (t))...\G N (t))+...+ CN (t)\G 1 (t))...\G N ^ 1 (t))\ s ), 

(10) 

where for simplicity we again assume kg = k w (no phase 
factors). The key point is that the coefficients c n (t) in Eq. 
(ITU1) are given by the single-atom calculation described in 
the previous section. Eq. © shows that c n (t) is in fact 
proportional to /3 n (t), which depends on the detuning of 
the corresponding (n-th) atom with respect to the laser. 
The collective atomic state is thus proportional to 

(3 1 (t)\ s )\G 2 (t))...\G N (t))+...+(3 N (t)\G 1 (t))...\G N ^ 1 (t))\ s ). 

(11) 

After the detection of the considered Stokes photon the 
states \G n (t)) will continue to evolve, and further photons 
will be emitted into other directional modes (let us recall 
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FIG. 3: We consider the situation where the energy of the 
excited state e is not the same for all the atoms in the en- 
semble, but has a Gaussian distribution with FWHM a. We 
have chosen a — 500 MHz in all our examples. We denote the 
detuning of the write laser with respect to the center of this 
distribution by Ao, which varies between and 1.25 GHz in 
our examples. 



that we are interested in the regime where the probabil- 
ity to emit another Stokes photon into the same mode is 
small). However, in the readout process the correspond- 
ing atomic excitations in s will just lead to the emission 
of additional Anti-Stokes photons in other, undetected 
directions. As long as the total number of excitations 
created in s is much smaller than the total number of 
atoms N, this has no significant effect on the collective 
interference in the readout process. 

Neglecting the additional excitations in s discussed in 
the previous paragraph, the starting state for the readout 
is given by 



Pi(ts)\s)\G 2 (tM))...\G N (t M )} + 
+PN{t s )\G 1 {t M ))...\G N - 1 {t M ))\s) 



(12) 



where ts is the time when the Stokes photon was de- 
tected, whereas £m is the memory time, i.e. the time 
when the excitation is read out, which may be much 
larger than ts- For long enough tu, depending on the 
lifetime of e, all the states \Gk) will essentially be equal 
to \g). The state of the atomic excitation in s is then 
proportional to 



P 1 (t s )\s)\g)...\g) 



■fa{t s )\g)..-\g)\t 



(13) 



We define the collectivity C, which is a function of the 
time of emission of the Stokes photon ts, as the fidelity 
of the (normalized) state of Eq. (fl"3"|) with respect to the 
ideal state of Eq. (fl|). 



IV. NUMERICAL RESULTS AND DISCUSSION 

We will now consider concrete examples. It is conve- 
nient to approximate the distribution of frequencies of 
the N atoms by a continuous spectral distribution n(A) 
with J e?An(A) = N. We will assume a Gaussian distri- 
bution with FWHM a for simplicity. We denote by Aq 



the detuning of the laser with respect to the center of the 
atomic distribution, see Fig. [3J 
The collectivity is then defined as 



G{ts) 



\fdAn(A)f3(t s ,A)\ 2 
N J dAn(A)\/3(tg , A) | 5 



(14) 



where ts is the time of emission of the Stokes photon 
as before. The relevant times are therefore those where 
Stokes photon emission is likely. As we have seen be- 
fore, the Stokes emission probability for an atom with 
detuning A at time t is proportional to the population 
in the excited state, \(3{t, A)| 2 . Let us recall that the 
total emission probability for the Stokes photon can be 
obtained by summing this quantity over all atoms (there 
is no collective interference for the Stokes photon emis- 
sion) . It is thus of interest to consider the average excited 
state population 



Pe(t) 



1 

N 



dAn(A)|/3(i,A)| 5 



(15) 



The values of ts that are likely to be observed are those 
for which p e (tg) is significant. 

For the inhomogeneous broadening we choose a = 0.5 
GHz, which is a realistic value both for hot gases 114 . 
EE EE EE HI and for rare-earth doped crystals [26j]. 
We furthermore choose Q(t) to have Gaussian temporal 
shape with a FWHM pulse duration of 0.1 [is, centered 
at t = 0.2/Lts in the figures, and with a maximum Rabi 
frequency of 1 MHz (corresponding to 2tt x 10 6 rad/s), 
a choice that is again motivated both by existing hot 
gas experiments and by potential future experiments on 
rare-earth doped solids. 

Figs. 2]and[5]showp e (t) and C(ts) respectively, for dif- 
ferent values of Ao. One sees the transition from the far- 
detuned (adiabatic) regime, where the population follows 
the Rabi pulse and the collectivity is high for emission 
times during the duration of the pulse, to the resonant 
regime, where the population no longer follows the pulse 
and the collectivity is low. The solid lines in the figures 
correspond to the case T = (negligible atomic decay), 
which is realistic for rare-earth doped solids, where typ- 
ical excited state lifetimes are of order 100 [is to 10 ms 
[26| . The dashed lines correspond to T = 5 MHz, which is 
a typical value for hot gases. One sees that the collectiv- 
ity is very similar in both cases, but it stays high a little 
longer after the pulse in the case with spontaneous decay. 
Intuitively, the decay of the collectivity after the pulse is 
due to the fact that Stokes emission for these times is 
dominated by non-adiabatic contributions (fluorescence) 
[l2T |. The spontaneous decay suppresses these contribu- 
tions compared to the adiabatic ones (Raman scattering) 
and thus enhances the collectivity. For very early times 
(ts < 0.02 [is) there are oscillations in the collectivity 
which we ascribe to a beating between fluorescence and 
Raman scattering. However, the Stokes photon emission 
probability for these times is exceedingly low, cf. Fig. 2] 

In the absence of other effects the collectivity of a given 
atomic excitation in s would remain unchanged once the 
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FIG. 4: The average excited state population p e as a function 
of time for different values of Ao. The probability of emitting 
a Stokes photon is proportional to this quantity. The detuning 
Ao of the laser from the center of the atomic distribution 
takes the values 0, 0.5 GHz, 0.75 GHz, 1.0 GHz and 1.25 
GHz from top to bottom. The solid lines correspond to a 
situation without spontaneous decay, the dashed lines to a 
decay rate F = 5 MHz. The write pulse is the same as in Fig. 
El The blue bottom-most line shows the profile of the square 
of the Rabi frequency, Q 2 (t). One can see that p e follows fi 2 
for large Ao, in analogy with the single-atom situation of Fig. 
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FIG. 5: Collectivity of the created atomic excitation in s as 
a function of the time of emission of the Stokes photon ts for 
different values of A , for T = (full lines) and T = 5 MHz 
(dashed lines). The pulse and inhomogeneous distribution 
are the same as in Fig. [4] The square of the Rabi frequency 
fi 2 (t) is shown as a fine blue line. One can see that a detuning 
Ao = 1.25 GHz is sufficient to have excellent collectivity for 
all relevant values of ts. 



excitation has been created via the emission of the Stokes 
photon. In practice it will decay on a timescale given by 
the spin coherence time characterizing the atomic ensem- 
ble, which can be very long for solid-state ensembles [2fj| . 
The spin transition g — s can also be inhomogeneously 
broadened in these systems, however the associated de- 
phasing can be compensated using spin echo techniques, 
which have already allowed the demonstration of coher- 
ence times as long as 30 seconds 



27| . In the case of hot 




gases the storage time is also limited by the motion of 
the atoms [II HOI]. 
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FIG. 6: The quantity n(A)|/3(f , A)| 2 at to = 0.2^s (at 
the center of the pulse) as a function of A, showing which 
atomic frequencies contribute significantly to the emission of 
the Stokes photon, and thus to the atomic excitation created 
in s. The laser frequency is defined to be A = 0. The center 
frequency of the atomic distribution takes the same values as 
in Figs. [4] and [5] The contribution from the resonant atoms 
(fluorescence) dominates for small values of Ao, but becomes 
negligible compared to the contribution from the bulk of the 
atomic distribution (Raman scattering) for large Ao. 



Fig. [5] shows the quantity n(A)\(3(t , A)| 2 as a func- 
tion of A for to = 0.2/is, i.e. at the center of the pulse. 
This permits one to see which frequency classes of atoms 
contribute significantly to the Stokes emission, and thus 
to the collective atomic excitation in s that is created, 
for the different values of Ao- One sees that for small Ao 
only the atoms that are resonant with the laser contribute 
significantly. On the other hand for large Ao the broad 
contribution from the bulk of the atomic distribution be- 
comes much more important than that of the resonant 
atoms (because there are so few of the latter). This is 
consistent with the results for the collectivity shown in 
Fig. [5] Intuitively, the sharp feature on resonance cor- 
responds to fluorescence, whereas the broad off-resonant 
contribution corresponds to Raman scattering [12j . 



V. CONCLUSIONS AND OUTLOOK 

We showed how one can quantify the collectivity of 
atomic excitations created by spontaneous Raman emis- 
sion in inhomogeneous ensembles, based essentially on 
simple calculations for a two-level system in combination 
with suitable averages over the inhomogeneous atomic 
spectral distribution. We found that quite moderate de- 
tuning (of the order of twice the inhomogeneous broad- 
ening) is already enough to be in the regime of high 
collectivity, where collective interference effects are as 
strong as in homogeneous systems. This is encourag- 
ing for future experiments in solid-state systems. It is 
worth mentioning that these results are in good corre- 
spondence with the conditions used in practice in hot 
gas experiments 0, [ID, [l|| ■ Of course the precise shape 
(not just the width) of the atomic distribution has to be 
taken into account for any given experiment. Note that 
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we have assumed that the inhomogeneous distribution is 
static, which is an excellent approximation for solid-state 
atomic ensembles. In hot gases collisions cause frequency 
changes, which become particularly important for longer 
write pulses. 

In this paper we have focused on the creation of an 
atomic excitation in s, i.e. the write process of the DLCZ 
protocol. The readout is a priori more complicated be- 
cause the Anti-Stokes emission exhibits collective inter- 
ference (and thus seems to be less amenable to a single- 
atom based treatment), and because the Anti-Stokes pho- 
ton can be reabsorbed by the ensemble, in contrast to 
the Stokes photon which couples to an essentially un- 
populated transition. However it should be possible to 
extend the present approach to a detailed study of the 
readout as well. The simplest case is a short, intense 
7r read pulse that excites all atoms from s to e simul- 
taneously, as described in the introduction. Once the 
excitation has thus been transferred to the excited state 
e, the problem is equivalent to the two-level situation 
studied e.g. in Refs. [H, HH and should thus be solv- 
able using the same techniques based on the Maxwell- 
Bloch equations for inhomogeneous systems. Note that 
efficient readout of excitations in e is possible in inhomo- 
geneous, absorbing systems even in the absence of con- 
trol beams, i f ap propriate phase matching conditions are 
fulfilled [13, |28|. These calculations are typically done 
in a one-dimensional approximation, which should well 
describe situations where write and read pulse, Stokes 



photon and Anti-Stokes photon all propagate along the 
same axis (in forward or backward direction). However, 
it may also be possible to extend the three-dimensional 
descriptions of Refs. [1(1 HI] to the inhomogeneous case. 

In practice, short tt pulses can be hard to implement 
due to laser power limitations and due to the risk of in- 
ducing unwanted transitions to nearby levels. It is there- 
fore of great interest to investigate alternative readout 
schemes using chirped pulses. This requires a detailed 
study of the impact of such excitation schemes on the 
phases in the collective excitation. EIT effects during 
the read pulse may also play a role in such a scenario 
[111 ]. EIT in inhomogeneous systems has been studied 
e.g. in Refs. [13, 

Let us finally note that inhomogeneous broadening can 
also have desirable effects. In the context of quantum 
memory protocols it should allow the efficient imple- 
mentation of temporal multiplexing [22]. Such "multi- 
mode memories" promise great speedups in the context 
of quantum repeater protocols [321 ]. It is a fascinating 
question whether a similar enhancement is possible for 
the DLCZ protocol. 
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